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E have recently measured the linear 

thermal expansion between room temper- 
ature and 950°C of ten samples of specially 
purified iron by means of a vacuum inter- 
ferometer.' The results for two representative 
samples (sample A, carbonyl iron held in 
hydrogen at 1500°C for 18 hours, and sample B, 
carbonyl iron melted in hydrogen) are given in 
Table I and Fig. 1. The temperature variation of 
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Fic. 1. Thermal expansion coefficient of iron as a function 
of temperature. 


the true expansion coefficient shown is typical of 
all but two of ten specimens examined; with 
rising temperature the coefficient increases, passes 
through a maximum between 400° and 500°, falls 
to a minimum at, or very near, the Curie point 


1 Austin, Physics 3, 240 (1932). 


then increases again to the a—y transformation, 
that is, to the change from a body-centered to a 
face-centered lattice. Quantitatively the coeff- 
cients for temperatures below 600° are in 
agreement but at higher temperatures there are 
differences, too great to be accounted for by 
experimental error, which indicate that between 
600° and 900° the thermal expansion is sensitive 
to very small differences in composition or to 
previous treatment or to both. 

There is no abrupt change in length at the 
Curie point, only a point of inflection on the 
length-temperature curve. No completely satis- 
factory explanation of the minimum coefficient 
has been found but calculations by Powell?* 
suggest that it is associated with the change in 
the strength of the binding between atoms which, 
on Heisenberg’s theory of ferromagnetism, 
accompanies demagnetization. 

The temperature of the transformation from 
body-centered to face-centered structure appears 
to be very sensitive to the presence of even slight 
amounts of impurity, varying with the sample 
from 903° to 940°. Sample A transformed at 
928+2° on heating and on cooling (rate ca. 
3.0°/min.), sample B at 915+2°. Judging by 
the behavior in passing through this change in 
lattice, sample A is the purest of all the specimens 
studied and we believe that 928°, although 20° 
higher than the value now widely accepted, is 
close to the true transformation temperature of 


2 Powell, Proc. Phys. Soc. London 42, 390 (1930). 
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TABLE I. True linear expansion coefficient a: =(1/l)(dl/dt) }. 
Sample 100° 200° 300° 400° 500° 600° 700° 800° 
A. Powdered carbonyl iron sintered in 
H; for 18 hr. at 1500°C 12.7 13.6 14.6 16.7 16.8 16.3 16.1 16.5 
B. Powdered carbonyl iron melted in 
Hy 12.7 13.8 14.6 16.6 16.8 15.8 15.2 14.9 
pure iron. In all but three of the specimens the eres ee ee ee ee 
change in length in passing from face-centered to mi ’ ] 
body-centered cubic lattice on cooling was Ey a J 
greater than the reverse change on heating. The ‘6h . 1 
cause of this discrepancy is not certain and is r 1 
being further investigated. "T | 
Ebert*® has shown that below 0° the thermal ak 4 
expansion of iron can be computed with reason- r » 
able accuracy from heat capacity data by means ‘rT | 
of Griineisen’s equation: eras een 
Q 100 200 300 400 500 600 700 800 900 1000 
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where 

Cur 
B= { "CT, a= ( ) 
0 (1/v)(dv/dT)7 ro 


and m and are exponents in the equation for the 
force between atoms. Taking Q)=170 kcal. 
(m+n-+3)/6=3.8 and deriving E and dE/dT by 
assuming a characteristic frequency of 400, all of 
which numerical values were used by Ebert, we 
find that the calculated coefficients for high 
temperatures are in very poor agreement with 
those observed, a result which is not surprising 
since it is well known that the heat capacity of 
iron at high temperatures cannot be represented 
by a single Debye function. The equation pre- 
dicts, however, that the ratio of a, to C, should be 





3 Ebert, Zeits. f. Physik 47, 712 (1928). 
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Fic. 2. Ratio of expansion coefficient to specific heat for 
iron plotted as a function of temperature. 


constant with temperature, a relation which is 
tested in Fig. 2 by using the expansion coefficients 
of sample A and the heat capacity data selected 
by Austin.‘ It is evident that proportionality 
continues with reasonable accuracy until the 
effects of the magnetic inversion begin to appear. 
Griineisen’s relation seems, therefore, not to be 
valid for ferromagnetics near the Curie point, at 
least insofar as macroscopic expansion is con- 
cerned. Whether true lattice expansion, deter- 
mined precisely by refined x-ray methods, would 
give closer agreement remains to be seen. 

The details of this investigation will appear in 
the Transactions of the American Society for 
Steel Treating. 


* Austin, Ind. Eng. Chem. 24, 1225 (1932). 
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An Inverse Boundary Value Problem in Electrodynamics 


L. B. SLICHTER, Massachusetts Institute of Technology 
(Received August 16, 1933) 


From knowledge of the electromagnetic field at the sur- 
face of a half-space due to a prescribed oscillatory source, 
is determined the unknown variation with depth of the 


conductivity and the dielectric constant. Unique solutions 
exist for these quantities. 





INTRODUCTION 


N THE inverse type of boundary value 

problem here considered, it is assumed that a 
given differential relationship is satisfied by a 
space field; that values of this field are prescribed 
over the boundary of a region; but that the 
differential relationships involve certain phys- 
ical properties of material which are unknown 
functions of position. The problem is inverse 
in that one seeks to determine unknown phys- 
ical properties of the region from a knowledge 
of the field on the boundary, rather than to 
obtain the field from a knowledge of physical 
properties. Such problems are fundamental in 
certain applications of physics, for example,' in 
the theory of geophysical prospecting, where it 
is desired to interpret the observed surface 
readings of a space field in terms of concealed 
geology. But, with a few exceptions, they are non- 
unique, and consequently the interpretation of 
the observational data is necessarily subject to 
basic uncertainties. However, unique solutions 
for two of the more simple problems of this 
character are available. Both, as might be 
expected, relate to a region in which the pertinent 
properties of the material are functions of one 
coordinate only, the depth. These two problems 
are (1) the interpretation, under the ray concept, 
of travel-time curves due to a seismic pulse; * and 


‘It appears possible that certain types of inverse 
boundary value problems will find application in the study 
of plasma. 

*See G. Herglotz, Phys. Zeits. 8, 145 (1907). Wiechert 
and Geiger, Phys. Zeits. 11, 294 (1910). H. Witte, Nach- 
richten der Gesellschaft der Wissenschaften zu Gottingen, 
p. 199 (1932). L. B. Slichter, Physics 3, 273 (1932). 


(2) the interpretation of a direct-current flow 
from a point electrode at the surface of a half- 
space of varying conductivity.* The present 
paper deals with the interpretation of the 
electromagnetic field produced by an oscillatory 
magnetic source situated at the surface of a plane 
earth whose conductivity and dielectric constant 
are unknown functions of depth.‘ 


§1.0. FORMULATION OF THE PROBLEM 


On the surface of a half-space whose unknown 
conductivity and dielectric constant vary only 
with depth, is located an electromagnetic source, 
consisting of concentric circular current sheets 
centering at the origin and varying in intensity 
with radial distance in a manner to be assigned. 
In time, the source currents vary harmonically 
with frequency w/27. The direct problem for a 
dipolar source and a constant conductivity is 
customarily set up in terms of a wave equation 
and a Hertzian potential function. However, it 
seems more convenient for present purposes to 
treat the problem in terms of the field equations 
directly, and for the case of a more general type 
of source. Cylindrical coordinates p, ¢, 2 will be 
used. 

It is clear from symmetry that the component 
B, of the magnetic vector is everywhere zero. 
Similarly, the components E, and E, of the 
electric vector must vanish at all points. Since 


3R. E. Langer, Bull. Am. Math. Soc. 39, 184 (1933). 
L. B. Slichter, Physics 4, 307 (1933). 

4 The corresponding direct problem for a dipolar source 
and an earth of constant conductivity has been treated by 
Sommerfeld and others. See A. Sommerfeld, Ann. d. 
Physik [4], 81, 1135 (1926). 
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the dependence of the vectors E and B upon 
time is sinusoidal, we set, as usual, 


E=R{\E, exp jut}, B=R{B, exp jot}, 

where E, and B, are now functions of position 
only. Then the two Maxwell circuital relations 
for the three nonvanishing components of the 
vectors EF, and B, are the following,® where, for 
convenience, these last subscripts will henceforth 


be omitted: 
4r 
-(- ser -(- -) == (oti) e, (1.1) 
Oz - we 4r 


OE 4/02 = +j(w/c)B, (1.2) 
8E,/ap +E./p =—j(w/c)B,. (1.3) 

In addition, 
div B=0; div «eE=4rp. (1.4) 


But the latter equations, (1.4), are automatically 
satisfied under the symmetry postulated when- 
ever the circuital relations are fulfilled. For 
simplicity, let it be assumed that »p=1. 

The boundary conditions to be associated with 
these equations are the following: At z=0, both 
FE, and B, are continuous; B, however possesses a 
discontinuity at z=0 determined by the surface 
currents of the source in accordance with the 
equation, 


(B,)+2—(B,)-2= (42/c)l(p), (1.5) 


where /(p) is the amplitude of the surface current. 
At infinity, the field vectors vanish at least as 
rapidly as in a spherically diverging electro- 
magnetic wave; that is, as 1/r’. 


$2.0. THe SOLUTION FOR Ey, 


After introducing the notations: 


5(z) =(o0(s)+je(s)w/4a)4a/c; g=jw/c (2.1) 
and substituting from Eqs. (1.2) and (1.3) into 


(1.1), we obtain 


*’ The customary Gaussian system of units will here be 
used. The conductivity, o and electric vector, E are in 
electrostatic units, the magnetic vector B in electromag- 
netic_units. 
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where E is now written for E, without confusion. 
The substitution E=R(p)Z(z) separates the 
variables and yields the Bessel’s equation, 


R(p)+(1/p)R(p)+(2—1/p?2)R(p) =0, (2.3) 
and the equation, 
Z—(— ew?/c?+jb(z))Z =0, (2.4) 
where 
b(z) =42wa(z) /c?. 
Above the ground’s surface, z<0, Eq. (2.4) 
takes the form. 
Z—(N—w*/c2)Z =0, (2.5) 
whence, if ko? =w*/c?, 
Z=exp {+(d—k,?)!2}. (2.6) 


Since E and B remain finite at p=0, the 
appropriate solution of Eq. (2.3) is R=J;(Ap). 
Let a solution of (2.4) which vanishes at infinity 
be denoted by Z;,(z, \). Further characteristics 
of this solution will be developed presently. For 
positive values of z, the solution for E will be 
denoted by E,, for negative z, by E_. Let us seek 
solutions of type, 


E,= f Ti(dp) Fi(A)Zi(z, Add (2.7) 
0 


-{ Ji(\p) Fo(d) exp [+(A2— ho?) 'z ]dd (2.8) 
0 


where the functions F,(A) and F2(A) are such as 
to secure the satisfaction of the boundary con- 
ditions and the convergence of the integrals. 
The continuity of E at z=0, independently of 
p, implies that F2(A) = F,(A)Z,(0, A), a relation 
which also assures the satisfaction of the con- 
tinuity condition on B,. The remaining condi- 
tion at the surface, namely Eq. (1.5), may be 
written 


(OE, 02)49— (0E_ 02)_0 = j+rul c’, 
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that is, 


- : 4rwl(p) 
{ Ji(Ap) Fi(A) {Z1(0, ¥) —Z1(0, A) (A? — Ro?) dA = 7 —— — (2.9) 
+9 


From this equation the required function F\(A) may always be determined, when the source function 
l(p) is given, by application of the Fourier-Bessel integral theorem. In the present case, we are 
particularly interested in a source which is of dipolar type at large distances. Within this limitation, 
then, we shall assign the function /(p) in accordance with analytical convenience. A suitable choice is 
suggested by the discontinuous integral,® 


af J (da) J,-1(Ap)dd = (p/a)’"; O0< p/a<l 
“é 


2; p/a=1 


0; p/a>1. 
Whence, if »=2, and if the source is such that 
l=p/a,0<p<a; l=}, p=a; 1=0, p>a; 
then the boundary condition (2.9) will be satisfied provided F,(A) is determined by the condition, 
jJo(da)(4rwa /c2) = F(A) {Z1(0, ¥) — Z1(0, A)? — he?) . (2.10) 
Thus, finally, the desired solution for E, is 


Anwa Jo(ra)Ji(Ap)Z1(z, A)dA 
E,=j ig Oo : (2.11) 
0 Z,(0, h) _ (A? — ko?)§Z, (0, d) 





C2 
Assuming for the moment the uniform convergence of the integrals, the following expressions for 


B, and B, are derived from Eqs. (1.2) and (1.3): 


4na 7” Je(ad)Ji(pr)Zi(z, A)dd 
Disies 0 eeomeaartntes (2.12) 
Cc 0 Z,(0, A) —Z,(0, A) (A? — ke?)! 


4a f J (ar) Jo(pr)Z1(z, A)Add 
0 


B,=-— ' 
c Z(0, ¥) —Z,(0, A) (A? — ky ?)! 


(2.13) 


As a check, this latter formula may be compared with Sommerfeld’s‘ result for a vertical magnetic 
dipole above an earth of constant conductivity. If o and € are constant, and k,? = eyw*/c* — j4rwo,/c’, 
then 

Zi=exp { —(—k,?))z} 
and Eq. (2.13) becomes 
4na 7” J2o(ar)Jo(pr) exp [ — (A? — k,?)!z Jada 
B,=— ~ -- ; 
Cc v0 (A? — ky?)! + (A? — hy?)! 


If Sommerfeld’s expression for the Hertzian potential is formally differentiated under the integral 
sign, a factor \? replaces the J,(ad\) above. Hence the integral diverges when z=0, and a check is 
clearly impossible. But for small values of \ and a, a formal resemblance between the formulae is 
apparent on account of the relation J2(ad)—a’d?/8. 


*N. Nielson, Handbuch der Theorie der Cylinder Funktionen, Teubner, 1904, p. 200. 
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§3.0 VALIDITY OF THE INTEGRALS FOR B, AND B, 


Granting, for the moment, the convergence of the formula (2.11) for E, it is desired to show that the 
integrals (2.12) and (2.13) do in fact represent the quantities B, and B, in accordance with their 
previous defining Eqs. (1.2) and (1.3); that is to say it is desired to establish the uniform convergence 
of the integrals (2.12) and (2.13). 

Firstly it will be shown that the integrands in Eqs. (2.12) and (2.13) are bounded for all p and z, 
and for all real, finite, positive values of X, i.e., that the quantities, 

rZ(z, r) | 1 —Z,(z, d) 1 
(A) ——-— a= --- — —— —_—— |, (B) : : 
Z,(0, d) |Z,(0, ¥)/Z(O, d) — (X22 — Ry?) Z,(0, d) (Z,(0, ¥)/Z1(0, ) — (2 — Re?) 


are bounded. 


To see this, set Z= R exp [1 fB(z, \)dz | in Eq. (2.4), with R, 8, real. Then 





R—(\?— ew?/2+6)R=0, (3.1) 
8+2(R/R)B—b(z) =0. (3.2) 


If, as it will be shown, the coefficient of R in Eq. (3.1) is always positive, then a positive mono- 
tonically decreasing solution of Eq. (3.1) is known to exist. The positive nature of this coefficient will 
first be postulated, and it will then be shown that a solution of Eq. (3.2) for 8 which insures the 
realization of this hypothesis can always be found. Then, since R>0O, Eq. (3.2) has no singularities 
except at infinity, hence possesses the general integral, 


1 “ 
B(z, d) -=—|s0)+ | 
9 d) Jo 


- o 
“9 


b(é) R°(E, nae (3.3) 


Now suppose for large z, €=«€2, ¢ =o2 (constants). Then 


Z\ =exp [ —2(A?— ke”)? ], (3.4) 
where 
ke? = €qw®/c? —j4 rwo2/c? 
or 
Zi = R2 exp 1825, (3.5) 
with 
R2=exp [ — { (A? — eqw?/c?)? + (4 woe/c*)?}* cos (y/2) ] (3.6) 
Bo= — | (A? — €2 w?/c*?)? + (42 woe/c?)*}! sin (Y/2). (3.7) 
Here 


tan Y =42wae/ (dc? — €ew"). 


As \ ranges from zero to infinity y may be chosen to decrease from ¥» <7 to zero. Thus sin (y, 2) and 
cos (¥/2) are always positive; 82 is always negative, vanishes as — 1, \ for large \, and is independent 
of s. Now the solution for 8(s, \), as expressed by Eq. (3.3) must approach 82(d) for large values of z. 
In the right member of Eq. (3.3), the numerator is an increasing function of z, since the integrand is 
positive. The denominator is a decreasing function, hence 8(s) steadily increased with z, and eventu- 
ally approaches the negative constant (as to 2), 82(A). Hence §(s, d) is always negative, and 6*(z, \) 
= B2*(A). Hence in Eq. (3.1) the coefficient \* — ew?/c?+8?=\ — ew, 'c? +82", where, from Eq. (3.7), the 
right member is a minimum for \=0. It follows that the coefficient is assuredly positive when e2 and 
o2 are so assigned that 
e(s) { e(s) —_ €2} <47°o2" ‘ow 


or, more simply, when ¢<es. Either condition may be considered to be satisfied at great depths 
without essential restriction on the physical problem. Thus both R and 8 decrease in absolute value 
as s increases; R remains positive, 8 negative, and B=». 
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Now in the quantity (B), the factor, 
Z(s)/Zs(0) = (R(@)/R(O) +i8(@)R@)/RO)} exp (i fad), 
0 
has an absolute value not greater than ((R(0), R(0))?+8°(0))!. Also, the absolute value of the 
denominator of the second factor, 
| Z1'(0)/Z1(0) — (A? —ke®)*| = | R(O)/R(O) +48 (0) — (2 ke?)}], 


‘ 


is at least as great as the absolute value of its imaginary component, | 8(0)| = |B2|. A similar argument 
obviously applies to the function (A). Hence 


‘vy —_ Ft 


- 


1S 
1e 





Ay 


A 
Jo(da)J;(Ap) im dx, i=0, 1 


0 
are proper integrals. 
When } is sufficiently large (say X=X,) then from Eq. (2.4) Z,;~e-* and the quantities (A) and (B) 
are asymptotically constants. Hence the integrals, 


[J (xa) Ji(A yf 4? ir 

a i ada, 
| te *"\ByI 
are essentially 


[ Jalar) Tu(oryean, ¢=0, 1. 
N 


1 


which converge uniformly for z=z,>0. Thus the integrals 


t (A) 
Ja{an) Hod) lan 
0 (B) 


have been shown to converge uniformly when z=2, >0, and the formula (2.12) is valid when z= 2, >0. 
Now 


2, (A) 
J2(anr)Ji( pr) ban 
J, \(B) 

is continuous at z=0, and it is known’ that the same is true of 


{ Jo(ar)Ji(pr)e**drd; pa. 
A 


Thus the right members of Eqs. (2.12) and (2.13) are continuous at z=0+-. Since the left members 
are also continuous at 0+, the equalities (2.12) and (2.13) have been established for 2-0. 

The similar integral (2.11) for E can be shown to converge by the same argument, so that the 
assumption introduced at the beginning of §3.0 is justified.* 


§4.0. INVERSION OF THE FORMULAE FOR B, AND B, 


From Eqs. (2.12) and (2.13), the expressions for B, and B, at the ground’s surface are 








B,(p, 0) = { Ji(Ap)k,(A)ddr (4.1) 
“0 
7G. N. Watson, Theory of Bessel Functions, Cambridge * The writer is indebted to Professor R. EF. Langer for his 


University Press, 1922, p. 399. helpful criticism of §3.0 in manuscript. 
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B,(p, 0) = f To(dp)k(A)dd, (4.2) 

where ° 
k,(A) =(4 ra/c)J2(da)Z3(0, r)/[Z1(0, dA) —Z1(0, A) (A? — Ro”)? ] (4.3) 
k() =(—4na/c)Jo(Aa)AZ1(0, d)/[Z1(0, d) —Z1(0, d) (2 — ho?) *] (4.4) 


are in general complex numbers which approach a real limit (see §3.0) as \ becomes large. Let the 
real and imaginary parts of the four quantities above be designated by notations of the type, 


R(X) =Rp' (A) +jk,'’(A), ete. 


If it be assumed that the quantities B,’, B,’’, B.’ and B,’’ at the surface are known by observation 
at all distances from the source, then the four functions k,’(A), -- -, R.’’(A), may each be determined by 
inverting Eqs. (4.1) and (4.2) in accordance with the Fourier-Bessel integral theorem.* This theorem 
states in effect that 


F(r) = f Juryudy { F(R)IuR)RAR, 
0 0 
provided 
i] F(R)R'dR 
0 


is absolutely convergent. 


Hence B,’, ---, B.’’ satisfy formulae of the type, 
By, o)= J “Tidoynan [~ BO, F109) de, (4.5) 
0 
provided the functions k,’, ---, k.’’ are expressed by the equations: 

k,' (A) =r f B,'(0, p)Ji(Ap) pdp (4.6a) 

0 
ke" O)=2 f By", oOo) ede (4.6b) 

0 
R/O)= [BL (0, p)Joldo)edp (4.6c) 

0 
ks"O)=r [- BL", 9)Joldo)edp (4.6d) 

0 


and provided that the integrals of type, 


f B,'(0, p)o!de, 
0 


are absolutely convergent. Since B is everywhere finite and vanishes at infinity as 1/p’, this latter 
condition is always satisfied. Thus it follows from Eqs. (4.3) and (4.4) that 


*G. N. Watson, Theory of Bessel Functions, §14.4, Cambridge University Press, 1922. 
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ke) Z(0, 0) fa By'(0, ») Jip) pd p+j So°B," (0, p)Ji(dp) pd 
kX) AZO») SoBz/(0, p)Jo(Ap) pdpt+j Sy?B.'"(0, p)Jo(dp) pdp 





(4.7) 


where the quantities occurring in the integrands may all be considered known. Yet it should be noted 
that the experimental observation of the real and imaginary parts of the B vector involves the 
difficult measurement of the phase of the vector with respect to the source current.'® However, it seems 
desirable to suppress questions of observational technique at present, since our primary object is the 
examination of the theoretical nature of the interpretation problem. We assume, then, that the 
function — Z,(0, \)/AZ(0, d) which we shall denote by K(A) and refer to as the “kernel,” is known." 
It is clear that this function embodies, and is precisely equivalent to, the totality of observational 
data at the boundary z=0. For, on the one hand, knowledge of B,(0, p) and B,(0, p) suffices to define 
K(aA). Conversely, knowledge of K(A) determines both B,(0, p) and B.(0, p) through Eqs. (4.1) to 
(4.4) inclusive. A convenient method due to Langer* may now be applied for determining the un- 
known conductivity function from knowledge of the kernel and the form of the differential Eq. (2.4). 





$5.0. DETERMINATION OF THE CONDUCTIVITYAND Also, in view of (5.1), (5.2) and (5.3), v(z, d) 


DIELECTRIC CONSTANT FROM THE KERNEL possesses an asymptotic expansion of the form, 
It is evident that the solution, Z;(z, 4), of Eq. oo 
(2.4), is of the form, v(z, A)~d 2 an(z)a", (5.5) 
Z1(2, A) ~e™?, (5.1) where 
ao= 1, (5.6) 


for large values of \. Consequently the kernel, 
and each coefficient satisfies the boundary con- 


K(\))=—Z,(0,d)/AZ(0,), (5.2) ition, 


a;(z)9a; as 2-0, (5.7) 
meee ee 9 ltt for large d and is analyticin the interval (0, ©). Substitution 
7 P » of the expression (5.5) for v(z, A) into Eq. (5.4) 
co leads to the relation, 
K(\)~ Dad", (5.3) 
n=0 co 
where do9=1, and the other constants are to be 2 Bak” = J(2) — bie), (8) 
considered known, since K(A) is known. They are 
in general complex. where the coefficients 8, are given by the formula, 
If the substitution v(z, \) = —Z,(z, A)/Z,(z, d) sil 
be introduced into the differential Eq. (2.4), Ba= —Gasit Dd. case-rOr (5.9) 
v(z, X) is seen to satisfy the Riccati equation, r=0 
b—v?+ {rn?—k,2+ jb(z) | =0, (5.4) But the right member of Eq. (5.8) is independent 


where of \. Hence each 8, except 8) must vanish. Thus 
ki? =c(z)w°c?. 


n+2 
P . Anpi= >» Qn+2—-rAry n+0. (5.10) 
10 Apparently the experimental determination of phase r= 
is most readily accomplished by transmitting the time sig- 
nals by an auxiliary high frequency wave. When the fre- For n= —1, this yields 
quency term is high relative to the conductivity term 
(see Eq. 2.4) the velocity of propagation is practically in- ay = 2a, =0 


dependent of the subsoil. 


11 The integrals occurring in Eq. (4.7) are of a type by Eq. (5.6), and consequently, for n= 0, we have 
readily evaluated by the Bush Mechanical Integraph; 


V. Bush, J. Frank. Inst. 212, 447 (1931). jb(z) —k:?(2) = 2ae(z). (5.11) 








418 ie B. 


Eq. (5.10) expresses the derivatives a;’(z) in 
terms of the functions a;(z) themselves. Indi- 
rectly it yields corresponding expressions for the 
higher derivatives, a; (z) as well. Now the 
functions a;(z) are all known at z=0, through 
the relations (5.7); hence the MacLaurin series 
expansion of each of the functions a; is com- 
putable. In particular, the first several coeffi- 
cients in the expansion for a2(z) are the following: 


a2(0) = de, 

c2(0) = 2a3, 
a2(0)/2!=2a4+a:2?, 
a2(0)/3!=4a5/3+8a3a2/3, 


ao (0) /4 ! = 2a¢/3 + 2a2a4+5a3?/3+2a2*/3, 


The function a(z) may thus be computed. By 
equating real and imaginary parts in Eq. (5.11) 
both the quantities b(z) and k,(z) are obtained; 
that is, the desired computation of the con- 
ductivity function as well as the dielectric 
constant has been effected. It is to be noted that 
the calculation of the conductivity function, 
b(z), from ae(z) necessitates definite knowledge of 
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the imaginary part of a2(z), and consequently of 
the real and imaginary parts, B’ and B”, in Eq. 
(4.7). In other words, if the observational data 
are restricted to points so near the source that 
the phase change is inappreciable, then, ipso 
facto, the computation of the variation of 
conductivity with depth is impossible. 

In the case of a symmetrically anisotropic 
medium, whose conductivity components o, and 
ao, vary only with 2, it is clear that only the 
horizontal component o, is effective in the 
present problem, since the secondary currents are 
confined to horizontal planes. 

It will be noticed that the interpretation prob- 
lem has been solved without explicit utilization of 
the frequency ; however knowledge of both B, and 
B, has been presumed. If, as is usual, the fre- 
quency w/27 is known, then either Eq. (4.3) or 
(4.4) suffices to define the kernel, 


—Z,(0, d)/dZ,(0, 2). 


It follows that B,(0) and BO) are mutu- 
ally dependent, and that knowledge of either 
alone serves to define the other and to furnish 
unique solutions for the variation of the con- 
ductivity and of the dielectric constant with 
depth. 
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A Study of the Evaporation of Small Water Drops 


H. G. HouGuton, Department of Electrical Engineering, Round Hill Research Division, 
Massachusetts Institute of Technology 


(Received August 21, 1933) 


A theoretical expression for the evaporation of small 
spherical water drops in still air is developed from the 
general evaporation equation of Jeffreys. Experimental 
data for the evaporation of drops ranging from 25 to 2600 
microns in diameter are obtained at several temperatures 
and relative humidities. After making approximate cor- 
rections for the cooling of the drops it is found that the 
results are in general agreement with the theoretical evapo- 
oration equation. Residual variations which are functions 


ys connection with the development of appa- 
ratus for the measurement of the liquid water 
content of a foggy atmosphere it was necessary to 
determine the time required for the evaporation 
of small water drops as a function of the drop 
size, the ambient temperature and the relative 
humidity. In view of the importance of this 
subject in other meteorological problems as well 
as in connection with air-conditioning and vari- 
ous industrial applications the results of this 
investigation are presented separately. 

This study was undertaken primarily for 
application to the evaporation of fog drops and 
since their velocity with respect to the air is 
small (0.1 to 2.0 cm/sec.) the problem may be 
considered as essentially equivalent to the evapo- 
ration of astationary drop into a still atmosphere. 
The case of larger drops falling freely in the air at 
a considerable velocity is more complex and is 
not considered in this paper. 

One of the first theoretical attacks on the 
general problem of evaporation was made by 
Stefan! who obtained analytical expressions for 
two special cases. One was for the evaporation 
from a tube of fixed length and constant cross 
section and the other was for the evaporation 
from a flush circular area, both into a still 
atmosphere. Much more recently Jeffreys’ has 


1Stefan, Sitzb. Akad. Wiss. Wien. 66, 385 (1873); 73, 
943 (1881). 
2H. Jeffreys, Phil. Mag. 35, 270 (1918). 


of the drop size and the difference between the water vapor 
density at the surface of the drop and at a distance from the 
drop are probably due to inaccuracies in the method used 
for computing the drop temperatures. The results of the 
paper furnish a means for computing the total time re- 
quired for the complete evaporation of a drop of liquid 
water into a still atmosphere at any given temperature 
and relative humidity. 


developed a more general solution which agrees 
essentially with Stefan’s results. A search of the 
literature has failed to discover either a theo- 
retical or experimental consideration of the 
evaporation of spherical drops but Jeffreys’ 
solution may be readily applied to this case. 
Jeffreys considers the evaporation problem as 
essentially one of gaseous diffusion. The general 
diffusion equation is 


dV @ aV 0 aV 0 OV 
LODO 26%). 0 
dt = ax Ox oy oy Oz Oz 


where V is the fraction of the air density at any 
point due to the water vapor (it is also called the 
concentration); ¢ is the time; d/dt is the total 
differential following a particle of the fluid; and 
k is the coefficient of diffusion. For the special 
case of evaporation into still air from a small 
surface it is shown that the diffusion equation for 
a steady state is given by 


vV?V=0, (2) 


where V?= the Laplacian operator. Jeffreys now 
observes that this equation is also satisfied by 
the potential in electrostatic problems and he 
proceeds to solve it by analogy to the electrostatic 
case. The solution for the total diffusion outward 


from the surface is found to be 
4nakCpo(Vo— V) g/sec., (3) 


where C is the electrostatic capacity of the 
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evaporating surface; p the density of the air; and 
V, and V the water vapor concentrations at the 
surface and at a considerable distance from the 
surface respectively. 

For the case of a sphere of diameter a the 
electrostatic capacity is a/2. If Vol is the volume 
of the sphere and 4 is the density of the water the 
evaporation equation becomes 


5-[d( Vol) /dt |= —2xkap(Vo— V) g/sec. (4) 


The rate of evaporation is thus proportional to 
the diameter of the sphere rather than to the 
surface. To facilitate computations it should be 
noted that the product of the concentration and 
the air density is the density of the water vapor 
and the above equation may be written 


5-[d( Vol) /dt]= —2xka(Do—D) g/sec., (5) 


where Dy is the water vapor density at the surface 
of the sphere; D the water vapor density at a 
considerable distance from the sphere. 

It was the purpose of the experimental work 
described below to verify this theoretical ex- 
pression and to indicate such modifications as are 
necessary. The drops were observed at rest 
suspended from fine wires or glass filaments 
within a small chamber in which the relative 
humidity could be maintained constant at any 
desired value. The entire apparatus was located 
in a temperature controlled room and the 
temperature changes were effected by varying the 
room temperature. The wires and glass filaments 
were made as small as possible in order to make 
their effect negligible. Since there is a definite 
lower limit to the size of drop which can be 
suspended from a support of a given diameter it 
was necessary to use supports of various sizes to 
cover the desired range of drop sizes. The di- 


ameter and material of the various supports and ~ 


the range of drop diameters with which they 
were used ‘are summarized in Table I. Over- 


TABLE I. 








Material of support 


Range of drop 
and diam. (microns) 


diam. (microns) 








Glass 250 over 2000 
Glass 85 2000 to 1000 
Glass 25 1200 to 300 
Platinum 

Wollaston 

Wire 2.5 500 to 25 
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lapping evaporation readings were obtained with 
various supports and it was found that the effect 
of the size of the support was inappreciable 
within the limits given in the above table. 
Drops much larger than 2000 microns were found 
to deviate from a spherical form and drops 
smaller than 25 microns could not be accurately 
observed because of the limitation imposed by 
the size of the smallest available support. Care 
was taken in all cases to insulate the supporting 
wires from the metal portions of the evaporation 
chamber so that the drops would be free to 
assume a temperature different from the ambient 
temperature. To reduce the area of contact 
between the drop and the support as much as 
possible the supports were given a very thin 
coating of paraffin. The drops were formed from 
distilled water by special pipettes with very fine 
outlets and were then transferred to the support. 
The smallest drops which could be handled 
satisfactorily were about 300 microns in diameter. 

The humidity within the chamber was main- 
tained constant by means of solutions of various 
inorganic salts which were contained in the 
bottom. In this way relative humidities of 0, 42, 
54, 66, 86.5, and 88.6 percent were obtained. A 
hair hygrometer was incorporated in the evapo- 
ration chamber to indicate when the humidity 
had reached a constant value and to detect small 
changes in relative humidity during a series of 
observations. The absolute values of the relative 
humidity maintained above the various solutions 
were determined in a separate chamber by means 
of an Assman psychrometer and a dew-point 
apparatus. A range of temperatures from 3° to 
32°C was used. 

The diameter of the drop under observation 
was measured by means of an ocular micrometer 
in a low power microscope. The drop was lighted 
from the side by a simple illuminator consisting 
of a 50 watt concentrated filament lamp, a single 
lens and a cell containing a ten percent solution of 
copper sulphate to remove the heat. The effec- 
tiveness of this cell was tested by placing a 
thermometer with a blackened bulb in the beam 
and over a period of one-half an hour the 
temperature rise did not exceed 0.1°C above the 
ambient temperature. As a further check on this 
important point two drops of equal size were 
allowed to evaporate for a given time, one with 
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the light on continuously and the other with the 
light off. The two drops were found to evaporate 
at the same rate within the limits of accuracy of 
the measurements. 

For a given temperature and relative humidity 
overlapping and duplicate data were obtained 
over a range of drop diameters from about 50 to 
2000 microns and in some cases the measure- 
ments were extended to include drops as small as 
25 and as large as 2600 microns. In general the 
procedure was to place as large a drop as possible 
on a given support and allow it to evaporate until 
it became so small that the support had an 
appreciable effect. Readings of the drop diameter 
were taken at intervals varying from one-half to 
ten minutes depending upon the rate of evapo- 
ration. 

A group of typical drop evaporation curves is 
given in Fig. 1. It will be observed that all of the 
curves are apparently parabolic. We may readily 
develop the theoretical equation of these curves. 
Since 

da/dt = (2/ma*)d( Vol) /dt, 


from Eq. (5): 


5(da/dt) = —(4k/a)(Do—D). (6) 
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Fic. 1. Experimental curves for the evaporation of small 
water drops into a still atmosphere. A, 20.8°C, 88.6 percent 
relative humidity; B, 19.6°C, 53 percent relative humidity; 
C, 21.7°C, 0 percent relative humidity; D, 4.8°C, 0 percent 
relative humidity; E, 20.3°C, 42 percent relative humidity; 
F, 22.0°C, 0 percent relative humidity. 





OF 


WATER DROPS 421 
Integrating and letting A be the initial drop 
diameter and putting 6=1 for water, 


a’ = A*—8k(D )—D)t, (7) 


which is the equation of a parabola. If 7 is the 
intercept on the time axis or, in other words, the 
time required for the complete evaporation of the 
drop, Eq. (7) may be written as 


a* = 8k(D,)— D)(T—2). (8) 


In Fig. 2 some of the data of Fig. 1 have been 
replotted on logarithmic paper with (7—?) as 
abscissas. The dashed line has a slope of exactly 
two and the experimental curves have very 
nearly the same slope indicating that they deviate 
only slightly from true parabolas in accordance 
with Eq. (8). This, in turn, shows that the drop 
evaporation is proportional to the diameter. 

It now remains to show that the evaporation is 
also proportional to (Do—D). This is most readily 
done with the aid of Eq. (6) by plotting a- (da/dt) 
against (D)»>—D). The values of da/dt were 
obtained directly from the experimental curves 
while Dy and D are respectively the water vapor 
densities corresponding to saturation at the 
ambient temperature and to the particular 
relative humidity and temperature of the evapo- 
ration chamber. These data are plotted in Fig. 3. 
The dots represent points obtained at a constant 
relative humidity (0 percent) and variable tem- 
perature and the circles are from the data taken 
at constant temperature (20°C) and variable 
relative humidity. The two sets of data do not lie 
on the same curve and do not lie on a straight line 
as would be indicated by Eq. (5). Furthermore, 
the slope is considerably less than 4k if k is given 
its usual value of about 0.25 cm?/sec. The 
explanation of these discrepancies seems to lie in 
the fact that the drop cools to a temperature 
considerably lower than the ambient temperature 
due to the evaporation. The equilibrium temper- 
ature of the drop is evidently determined by a 
balance between the heat loss due to evaporation 
and the heat gained by radiation, conduction and 
convection. The exact analytic solution of this 
problem is, however, practically impossible be- 
cause of the multiplicity of factors which are 
involved. The actual measurement of the temper- 
ature of drops as small as were used in this 
investigation is also difficult and was not at- 
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Fic. 2. Some of the data of Fig. 1 replotted on logarithmic paper with (7-?) 
as abscissas. A, 21.7°C, 0 percent relative humidity; B, 20.3°C, 42 percent 
relative humidity; C, 19.6°C, 53 percent relative humidity; D, 20.8°C, 88.6 


percent relative humidity. 


tempted. A very similar problem which has 
received considerable attention is the measure- 
ment of the relative humidity by means of the 
temperature depression of a wet bulb ther- 
mometer. The so-called psychrometric formula 
has been developed for this case on the basis of 
certain theoretical considerations but it contains 
sO many approximations and experimentally 
determined constants that it is little better than 
an empirical expression. However, considerable 
experimental work has been done on this subject 
and the results for ventilated wet bulb ther- 
mometers are probably quite accurate. Data for 
unventilated wet bulb thermometers are not in 
good agreement possibly because of differences in 
the experimental conditions since the thermome- 
ter would seldom be entirely unventilated unless 
enclosed in a chamber. In view of these vari- 
ations, a wet bulb thermometer was prepared and 
inserted in the evaporation chamber in the 
position normally occupied by the drop. The 
temperature and relative humidity were then 
varied and the depression of the wet bulb 
thermometer noted. From these data the constant 
in the psychrometric equation was determined. 
For all values of relative humidity and ambient 
temperatures from 5° to 25°C the following 
equation holds quite well as long as the tempera- 
ture of the wet bulb does not drop below 0°C: 


e—e’=12.3X10-'*B(r—7’), (9) 


where 7 is the ambient temperature in °C; 
r’, the temperature of the wet bulb thermometer 
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Fig. 3. The relation between the apparent density differ- 
ence and the product of drop diameter and rate of change 
of drop diameter. 
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in °C; e, the vapor pressure in the chamber; 
e’, the saturation vapor pressure corresponding 
to r’; and B, the barometric pressure measured in 
same units as e and e’. For temperatures higher 
than 25°C the actual depression of the wet bulb 
thermometer was found to be somewhat greater 
than that given by the above expression. Al- 
though it is evident that these results cannot be 
expected to apply strictly to the water drops 
which are much smaller and of a different shape 
than the thermometer bulb and which are sus- 
pended in a different manner they will at least 
permit approximate drop temperatures to be 
computed. 

Corrected drop temperatures were computed 
for all of the drop evaporation data on the basis 
of this psychrometric equation. Dy then becomes 
the density corresponding to saturation at the 
corrected drop temperature. The data of Fig. 3 
have been corrected in this manner and replotted 
in Fig. 4. It will be observed that all of the points 
lie very nearly on a straight line. The slope of this 
line is 0.72 X 108u?/sec. giving a value of k=0.18 
cm?/sec. This is somewhat lower than the 
accepted value of the diffusion coefficient. As 
given in the International Critical Tables, k is 
about 0.22 cm?/sec. at 0°C and 0.25 cm?/sec. at 
20°C. Possible causes of this deviation are con- 
sidered below. 

Residual variations with a definite trend still 
remain after the correction for drop temperature 
has been made. Values of k have been computed 
from the data and plotted against drop diameter 
for various values of (D)—D) in Fig. 5. It will be 
observed that the apparent value of k decreases 
for small drop diameters and for iarge values of 
(Do—D) and that the maximum value corre- 
sponds very closely to the value of the diffusion 
coefficient at 20°C. This would seem to indicate 
that the actual drop temperatures were lower for 
small drops and large density differences than the 
temperature of the wet bulb thermometer. In 
view of the obvious differences in the two cases 
this discrepancy is not surprising. It might be 
possible to obtain a more satisfactory empirical 
expression for the drop temperature from the 
data of Fig. 5 by assuming a fixed value of k but 
this has not been done because it is not certain 
that the observed variations are entirely due to 
inaccuracies in computing the drop temperature. 
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Fic. 4. The data of Fig. 3 replotted against a density differ- 
ence obtained from the corrected drop temperatures. 
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Fic. 5. The apparent diffusion coefficient as a function 
of drop size and density difference. ((Do—D) values are 
x 10-6.) 


If 7 is the total time in seconds which is 
required for the complete evaporation of a drop of 
initial diameter A it follows from Eq. (7) that 


T = A?/8k(Dy—D). (10) 


To compute from this expression it is first 
necessary to determine the corrected drop tem- 
perature from Eq. (8). Dy and D may then be 
obtained from tables (density compared to 
water). An average value of k is then obtained 
from Fig. 5 and if A is expressed in centimeters 
Eq. (9) will give the total time in seconds 
required for the complete evaporation of the 
drop. 
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It has been noted earlier in this paper that this 
work was undertaken primarily for application to 
the problem of the evaporation of fog drops. The 
smallest drops which have been studied were 
about 25 microns in diameter and many fog 
drops are smaller than this.’ It does not seem 
practicable to extend drop evaporation measure- 
ments to smaller drops by the methods of this 
paper but it is hoped that confirmatory data for 
small drops can be obtained in the near future by 
another method. In the meantime it is not 
thought that serious errors will be introduced if 
the curves of Fig. 5 are extended to include drops 
having a diameter as small as 5 or 10 microns. It 
is probable that the larger fog drops (20 to 40 
microns) will evaporate somewhat more rapidly 
than the results of this paper would indicate 
because of their velocity with respect to the air. 
In connection with the application of the results 
of this paper to the evaporation of fog drops it is 
also important to remember that these measure- 
ments have been made on drops composed of 
distilled water whereas fog drops are usually 
composed of dilute salt solutions. For the fogs 


+H. G. Houghton, Physics 2, 467 (1932). 
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which are encountered along the southern coast 
of New England this salt content varies from 20 
to 400 parts per million by weight. This salt 
content will have little effect at first but as the 
drop continues to evaporate the concentration of 
the salt will increase and will reduce the satu- 
ration vapor pressure and so retard the evapo- 
ration. If the relative humidity is sufficiently high 
a point will finally be reached where the vapor 
pressure at the surface of the drop will equal the 
vapor pressure in the air and evaporation will 
cease. In most cases this effect will probably be of 
little importance since the size of the drop when 
equilibrium is reached will usually be not much 
greater than the irreducible minimum size repre- 
sented by the solid salt and dust which remain 
after complete evaporation. An effect which is 
very similar to the above will occur with very 
small drops of even pure water. This is due to the 
lowering of the saturation vapor pressure as the 
radius of curvature becomes very small. This 
effect is entirely inappreciable for drops greater 
than one micron in diameter and does not become 
of real importance until the drop is about one- 
tenth of a micron in diameter. 
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The Final Distribution of Retained Liquid in an Ideal Uniform Soil 


W. O. Situ, Pittsburgh, Pennsylvania 
(Received August 24, 1933) 


If liquid, initially saturating the pores of a soil, is al- 
lowed to drain from that soil, no such complete drainage as 
occurs in a simple capillary tube is observed; large quanti- 
ties of liquid are retained. The quantity of liquid held by 
the soil immediately after drainage gradually alters; it de- 
creases to a stationary value and the liquid so held finally 
arranges itself in a stable distribution. Whatever be the 
distribution, the bounding liquid surfaces are capillary 
ones and subject to the laws of thermodynamics and 
capillarity; the Kelvin relation, in general determines the 
curvatures. The problem is considered for an ideal soil, 
that is an assemblage of spheres, of a single size, packed at 
random. There are in the ideal soil three types of distribu- 
tion. First, there is a pendular region where the liquid is 
retained in the form of single rings of liquid wrapped 
around the axis of each pair of adjacent grains in contact. 
Second, there is a funicular region, which may be con- 
sidered to arise by coalescence of the rings; we have two 
or more grain contacts imbedded in a single liquid mass, 
that is, webs of liquid enmeshing two or more grain con- 
tacts and which may involve many grains of the packing. 


There is, finally, a saturation region in which all pores are 
completely filled. The limits of each of these zones are de- 
termined. The funicular zone is a zone of hysteresis and all 
values between a normal minimum and complete satura- 
tion are usually possible. To calculate the normal minimum 
distribution the actual packing is replaced by an average 
packing of grains in hexagonal array and equally spaced 
to give the required porosity. The funicular distribution is 
now in the form of single rings wrapped around the axis of 
each pair of grains of the packing. The volume of each type 
of ring is determined in terms of (o¢/pghr), where o is the 
surface tension of the liquid and p its density, r the grain 
radius, and h the height of the ring above a free liquid 
reference plane. The average volume of each type of ring is 
obtained for all rings confined between two given planes. 
From the average ring volume and the number of such in 
each cc of packed space, the total retention in each zone 
above saturation is found. The distribution, that is, the 
liquid retained in any lamina at a height h above the free 
liquid, is also given. The results are compared with the data 
of King. 





T IS well known that a long column of soil, if 
its pores have been initially filled with liquid 
and then allowed to drain, retains considerable 
quantities of liquid; and no such complete 
drainage as occurs in a simple capillary tube, 
subjected to a similar process, takes place. If the 
soil column has been rapidly drained, large 
quantities of liquid, distributed throughout its 
entire length, will be found. If, now, the rapidly 
drained column is allowed to stand for several 
months, and then examined, the distribution of 
liquid will be observed to have undergone a 
great change; the upper parts of the column will 
be found to contain but the slightest traces of 
moisture, while the lower parts still retain con- 
siderable amounts. If the column is allowed to 
stand further, and examined at intervals, the 
distribution will be observed to become finally 
stationary. We may, on the other hand, suppose 
the drainage to have occurred very slowly and to 
have occupied a time equal to that covered by the 


whole of the experiment just described; the 
distribution will again become stationary, but 
the total retention may be less than that for 
rapid drainage. In either case, the distribution 
finally prevailing depends solely on the vapor 
pressures of the liquid, existing in the pore space 
during the final stages of the process; it is 
permanent and thermodynamically stable. We 
shall consider the final distribution for an 
assemblage of spheres of a single size, packed at 
random, that is, an ideal soil. 

We shall, in order better to understand the 
nature of the final distribution, consider that the 
drainage of the ideal soil is performed under 
definitely controlled conditions. We suppose the 
soil to be packed in the column shown in Fig. 
la, in which it is supported by the gauze c. The 
system is kept at constant temperature and 
contains, besides the soil, only the liquid and its 
vapor. We suppose, finally, that all pores have 
been filled and drained, and that the system, 
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Fic. 1a. Single drainage column: J marks the free liquid 
surface; ¢ is gauze supporting the soil; Amin. marks the 
minimum height of rise and lower limit of saturation; hp 
denotes the lower limit of the pendular zone. 

Fic. 1b. Double column to show saturation hysteresis: 
hmin. Marks the minimum height of rise and Amex, the 
maximum rise; the intervening interval is the region in 
which saturation hysteresis occurs. 

Fic. 1c. Drainage column of King: B is the soil; D the 
collecting bottle; C and E are small capillaries. 


with the free surface of the liquid in a position 
just below the soil—marked J in the figure—has 
been allowed to come to equilibrium. It is 
convenient, now, to explore the drained column, 
throughout its length, and ascertain just how the 
liquid is distributed; whether in discrete masses, 
a continuous body or perhaps both. We shall 
consider this question in terms of the governing 
physics; the conclusions reached are easily con- 
firmed by experiment. 

If we examine the liquid retained in the upper 
parts of the soil column, we find it is in discrete 
masses and each mass is a ring of liquid wrapped 
around the contact point of a pair of adjacent 
grains as in Fig. 2a. Liquid must accumulate in 
these positions rather than in a uniform layer 
over the grain surfaces, otherwise the energy of 
the retained-liquid system would not be a 
minimum; the wetting layer, one molecule or so 
thick, is, of course, supposed to exist. As we 
proceed downward the rings increase in size until, 
finally, they begin to coalesce and merge into 
more complicated masses; we have, then, two or 
more contact points imbedded in a common 
liquid mass. 

The single rings do not merge at a definite level 
of the column, but do so progressively as we 
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Fic. 2. (a) Simple pendular ring of liquid. (b) Simple ring 
of liquid in spaced hexagonal packing. 


descend further down. That this is so may be 
seen from the fact that the larger the angle 
between the axes of two adjacent pairs of grains, 
the larger the rings that can be supported; the 
axis of a pair of grains is defined as their line of 
centers, and is coincident with the axis of the 
‘ring. Now the angle between the axes of two 
adjacent grain pairs varies within definite limits 
because of the random nature of the packing; 
thus coalescence of rings begins first with those 
grain pairs whose axes are closest, and occurs 
lower in the column with rings whose axes make 
the larger angles. The more complicated masses 
grow as we descend until finally a stage is 
reached where they merge and form a capillary 
surface extending completely across the soil 
body ; the pores below this surface are completely 
saturated. 

We thus observe, as we descend the soil 
column, three types of liquid distribution. First, 
there is a stage with single liquid rings wrapped 
around grain pairs, with one contact point to a 
single mass—called by Versluys! the pendular 
stage. Second, there is a stage where the single 
rings are combined, in general, into more 
complicated forms, and characterized by having 
two or more contacts imbedded in a single liquid 
mass—Versluys! calls this the funicular stage. 
Finally there is the stage of complete saturation 
in which a single mass of liquid, completely 
filling all pores, extends across the soil body and 


1 J. Versluys, Die Kapillaritdt der Boden, Inst. Mitt. f. 
Bodenk. 7, 117-140 (1917). 
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possesses a simple continuous capillary surface as 
boundary—it is the saturation stage. 

It is fundamental to observe that whatever the 
boundary surfaces of the liquid found in the soil 
body, they are capillary surfaces possessing 
definite curvatures and as such are subject to the 
laws of thermodynamics and capillarity. When 
the system is in stable thermodynamic equilib- 
rium, the vapor pressures have definite values 
throughout the soil column. In the experiment 
just described the vapor pressures at any two 
points, at an equal distance above the free liquid, 
are identical. Kelvin’s? investigations furnish a 
criterion for determining the curvatures of a 
stable capillary surface over which a given vapor 
pressure exists. 

Kelvin’s theorem enables us, further, to assign 
definite curvatures to surfaces in equilibrium at a 
given distance above a free liquid. For, according 
to Kelvin, if two similar masses of liquid, 
confined by identical boundaries, are at equal 
distances above a free liquid, and if, further, one 
is in communication with and the other isolated 
from the free liquid, then not only the vapor 
pressures over each, but the respective curvatures 
of corresponding parts of each surface are equal. 
We may, for example, consider a single ring of 
liquid wrapped around the contact point of a 
pair of grains and isolated from the free liquid, 
and a similar ring in communication with the 
free liquid—communication can be accomplished 
by a pipe, one end of which opens into the free 
liquid and the other, running through one of the 
grains, opens into the ring at the contact point; if 
both rings are at an equal distance h above the 
free liquid, the vapor pressures are equal and 
hence the curvatures are also equal. We can, 
therefore, consider a single isolated liquid mass, 
for example the single liquid ring, as if it were in 
communication with a free liquid and determine 
its equilibrium from the well known curvature 
relation of Laplace; the curvatures of all capillary 
surfaces existing in the soil body can be found. 
We shall, throughout this discussion, suppose 
complete wetting. 

The curvatures of any element of the capillary 


surface are given by the equation, 
Ap=o(1/R:+1/R:), (1) 


2 W. Thomson, Phil. Mag. [4] 42, 448 (1871). 
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where &; and 2 are the principal capillary radii 
of curvature of the given element, Ap the pressure 
drop across it and o the surface tension. If we 
multiply Eq. (1) by dS cos ¢ and integrate over 
the whole capillary surface, we find’ 


[Ap dS cos = fo cos 6 ds (2) 
where dS is an element of area of the meniscus, ¢ 
the angle which the normal to this element makes 
with the vertical z-axis, ds an element in the 
meniscus grain contour and @ the angle between 
the vertical and that normal to ds which is 
tangent to the meniscus. We shall apply Eq. (2) 
to a single ring of liquid, and then calculate the 
quantity of liquid retained in the pendular region 
of the soil. 

If we suppose complete wetting and apply Eq. 
(2) to the ring shown in Fig. 2a, we find, neglect- 
ing small quantities of the order of a grain radius 
that 


Ap2rr sin 0(2r —2r cos 0) = 2(2zr sin @)o cos 8, 


from which 


sec 0=1+0/(rAp), (3) 


r being the grain radius. The angle @ may be taken 
as a measure of the size of the ring; the angle of 
the ring is defined as the angle between the axis 
(line of centers) of the two adjacent grains about 
which it is formed and the radius vector from one 
of the grain centers to the outer edge of the ring 
on that grain, and is easily seen to equal @. Eq. (3) 
enables us, once the grain radius 7, the surface 
tension o¢ and pressure drop Ap are known, to 
fix definitely the size of the ring. For the purpose 
of the present problem we may set Ap=pgh, 
where p is the density of the liquid and g the 
acceleration of gravity; hence, we find 


sec 0=1+0/(pghr), (4) 


for the size of the stable ring at a height h above 
the free liquid. This equation is applicable so long 
as the single rings do not overlap and merge. 
To find the actual volume of the ring around 
two grains in contact we shall make the as- 
sumption that the section of the capillary surface 


3 See any treatise on capillarity; e.g., Webber and Gans, 
Repertorium der Physik. 
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made by a plane passing through the axis is the 
arc of the circle of radius R tangent to both 
grains, and that the volume is generated by 
revolving this arc around the axis. We shall let y 
be the radius of the circle traced by the capillary 
surface in a plane normal to the axis and midway 
between the two grain centers (cf. Fig. 2a). The 
volume v of liquid in the ring is not difficult to 
find; for example, let w, be the volume of the 
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figure of revolution formed by revolving the 
circular arc of radius R around the axis, and 
included between two planes normal to the axis 
and passing respectively through the points, on 
each grain, of tangency to the generating arc, and 
let we be the volume of the grain portions 
included between these planes, then v=w,—w». 
After a simple integration,‘ we find, on setting 
a=R/r and k=y/r,® 


v/(2mr*) =[a/(1+a) ][(k+a)*—a(k+a)}1—1/(1+a)*}! 
—a(k+a)(1+a) sin“ {1/(1+a)}+a?—a/(1+a)]. (5) 


From the geometry of the ring shown in Fig. 2a, we have the relation (1+ a)? = (k+a)*?+1, which may 


be used to simplify Eq. (5) to 


v/(2mr*’) =o?[1—(k+a) sin {1/(1+a)}]. (6) 


We may express & and a in terms of 6; from Fig. 2a we see that 


(1+a) = sec 0; 


(k+a)= tan 6; 


sin“ [1/(1+a)} =(2/2—8). (7) 


Substituting these expressions in Eq. (6) we find as an alternative form: 


v/(2nrr*) = (sec @— 1)? 1—(2/2—6) tan 6]. (8) 


Comparing Eqs. (4) and (7) we see that 
a=a/(pghr); 


k=[20/(oghr) + (o/pghr)? }!— (o/pghr). (9) 


If, finally, these values of k and @ are substituted in Eq. (6), we find 


v/(2mr*) = (o/eghr)*[ 1 — {20/eghr + (o/oghr)*}* sin {1/(1-+0/oghr)} J. 


One obtains from Eq. (10) the volume of liquid of 
density p and surface tension o in the form of a 
ring around the contact point of two grains each 
of radius 7 at a height / above a free liquid, and in 
stable equilibrium; g is the acceleration of 
gravity. Eqs. (4) and (8) may also be used to 
derive Eq. (10). Fig. 3 shows the quantity 
v/(2xr*) represented as a function of h. 

In order to compute the total retention in a 
soil column included between two planes at 


4 We find: 
w,/2e = {(y+R)*+R*} Roos 0—(y+R)R*(4/2—8) 
— R*(y+R) sin 6 cos 0—(R cos 6)3/3 
We/2n = (2/3)r+Pr (cos 6)?/3—r cos 8, 


where sin 6=(R+y)/(r+R) and cos @=r/(r+R). 
5 cf. Smith, Foote and Busang, Phys. Rev. 36, 524-530 
(1930). 





(10) 





distances h, and hy respectively above a free 
liquid, it is convenient to have an expression for 
the average ring volume (@) for all rings between 
the two heights A, and he; (h2>h,). We require 


v 1 he ¥V 
= f adh. 
2nrr*® he—hydn, 2xr' 


For this purpose we use Eqs. (4) and (8); we find 
from Eq. (4): 








(11) 


(—o/pgr)(sec @ tan @) 
(sec 6—1)*dé 





With this value of dh, we find from Eqs. (8) and 
(11): 





(12) 


B 1 —@¢ 
= f — ¥(6)dé, 
2nr® he—hy/Yr = pgr 
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Fic. 3. The quantity (v/27r*) represented as a function 
of h for Eq. (10). o/pgr=1. 


where 
¥(0) = {1—(w/2—8@) tan 6} sec @ tan 0. 
We set F(h) = {y(0)d6, and obtain® 


F(h) =} (sec 0) {1—(a2/2—6) tan 6} 
+3(a/2—6) log tan x+x{ (log x) —1} 


+ (x3/9){1+7x2/50+62x!/2205+---}, (13) 


where x=}3(27/2+0). To find F(t) and Fhe), 
values of 6, and 6 are determined by use of Eq. 
(4) from the corresponding values of fh; and In. 
We find from Eqs. (12) and (13) that the average 
volume of the ring in the pendular zone between 
the planes at distances h, and h, from a free liquid 
surface is 


6 We find 


SW(0) do— }(sec 6) {1—(2x/2—6) tan 0} 
+3S (x/2—6) sec 6 dé. 
Now 
JS (x/2—6) sec 6 dé =(x/2—6) log tan 3(x/2+80) 
+S log tan}(4/2+6) dé. 


To find Sf log tan }(1/2+6) d@ we set x= }(x/2+6) and 
integrate in series. The series for log tan x (cf. B. O. Pierce, 
Short Table of Integrals) is convergent for all values of 
x? <(x/2)?; @ is in general 30° or less, as will be shown later; 
hence x= }(x/2+6) generally has its greatest value near 
unity, and the integration is within the region of con- 
vergence. 
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D _(~a/pgr){ Fhe) — F(a) } 
" he—hy , 





(14) 
2nr* 


The limits of the pendular zone must be con- 
sidered. This zone is confined to the part of the 
soil column in which single rings of liquid exist; 
when the single rings begin to merge we pass into 
the funicular zone. It is not difficult to fix the 
value of @ at which the rings commence to merge. 
The soil, here considered, is composed of spherical 
grains of one size packed at random. The value of 
6 at which two rings touch and merge is deter- 
mined by the angle 6 between the axes of the two 
adjacent grain pairs supporting them. These 
relations are shown in Figs. 4, a, b, c. Obviously 
this phenomenon will first occur when the axes of 
the three grain pairs involved form an equilateral 
triangle; hence @ = 30°. We find from Eq. (4) that, 
for this value of 6, overlapping commences when 
(pghr/c) =6.452. Thus h=6.452(¢/pgr) fixes the 
lower limit of the pendular zone; the upper limit 
is the top of the soil column, including h infinite 
as the extreme upper limit. Values of (pghr/c) 
<6.452 down to values of this quantity required 
for saturation are in the funicular zone. 

To find Q, the total volume of liquid per cc of 
packed space, in the pendular zone, it is necessary 
to know the total number of contacts N per cc (of 
packed space). For this purpose we suppose the 
quantity N plotted, as in Fig. 5, as a function of 
the corresponding porosity P. N, and N, are 
points of this curve, where N, and N, are the 
numbers of contacts per unit volume for close 
hexagonal (P=0.260) and cubic (P=0.476) 
pilings respectively ; N;, = 6(2!)/8r’ and N,=3/8r', 
for two grains are required for each contact. We 
make, as a first approximation, a simple inter- 
polation between these values to find N. From 
Fig. 5 


(N-N.)  (Na-N) 
(0.476—P) (0.476—0.260) 





We find 
N=N,x+N(1—x), 
where 


x =(0.476—P)/(0.217), (15) 


and substituting the values of N, and N, given 
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Fic. 4. Transformation of single liquid rings. (a) Single rings with angle 5>60°. (b) Single rings about to 
merge in three grain pore with 5 =60°. (c) Simple web of liquid existing after coalescence of rings. (d) Rings 


of liquid in spaced hexagonal packing, with 6 = 60°. 


above we find’ 


N =3(1+1.828x) /8r°. (16) 


In the pendular zone, finally, the average 
volume of liquid Q, retained per cc of packed 
space is given by 


2,=iN, (17) 


where N is given by Eq. (16) and 6 by Eq. (14). 
In a volume V of packed space the total retained 
liquid is 2,V ce. 

For the purpose of determining the distribution 
of liquid in the soil, the quantity retained per cc 
of packing is required to be known for a lamina of 
soil at a height # above the free liquid. The 
variation of v in a small interval is approximately 
linear and we may, for this purpose, disregard the 
averaging process. We may set 


2,=vN, (18) 


where v is given by Eq. (10) and JN, as before, by 


7It is not difficult to show that the average number of 
contacts per grain derived from Eq. (16) is identical with 
the value given in an earlier paper (Smith, Foote and Bu- 
sang, Phys. Rev. 34, 1271-1274 (1929)). We note that the 
number of grains m per unit volume, in a packing of 
porosity P, is found by a similar interpolation. We find, as 
above, m=m,x+m,(1—x) where m, and m, are the num- 
ber of grains per unit volume for close hexagonal and 
simple cubic pilings respectively; m,=(24)/8r> and 
m,.=1/8r°. We find 


m =(1+0.414x) /8r. 


The average number of contacts per grain is given by 
n=2N/m; hence 


n = 6(1+1.828x) /(1+0.414x), 


which is the expression given in the earlier paper and found 
to agree with experiment. 
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Fic. 5. Number of contacts per cc of packing represented as 
a function of porosity P. 


Eq. (16). The value of v at the center of the 
lamina under investigation only need be con- 
sidered; the increased moisture below this point 
is closely compensated by the decrease in 
moisture above. 

The distribution of liquid in the funicular zone 
is complicated and presents, throughout, a 
variability not found in the simpler pendular 
zone; hysteresis is a commonly observed effect. 
Not only single rings of liquid but more compli- 
cated forms, namely, webs of liquid enmeshing 
two or more grains, are present. It is possible to 
have single pores completely saturated or even 
several associated pores in that state. For 
example, the random tetrahedron formed by the 
centers of any four associated grains, not in 
regular array, may be completely filled; or else, 
the retained liquid may be confined to masses 
symmetrical about the faces of the tetrahedron 
(i.e., the webs of liquid which may extend to 
adjacent tetrahedrons) ; or we may have, finally, 
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simple rings about contact points. The manner of 
drainage is usually the controlling factor. Before 
discussing these effects further, it is convenient to 
fix the limits of the funicular zone in which they 
commonly occur. 

The upper limit of the funicular zone is its 
junction with the pendular zone and, as stated 
earlier, is given by the value of h for which 
(pghr/o) = 6.452. 

The lower limit is given by the least value of h 
at which saturation occurs, that is, the minimum 
height of capillary rise. Liquid, rising into a soil, 
would ordinarily come to equilibrium when the 
meniscus occupies the widest opening available in 
its pore cells. Condensation of vapor, just ahead 
of the advancing meniscus, gives rise to funicular 
masses of liquid which pull the meniscus further 
into the soil; equilibrium is reached, when the 
rising liquid is just unable to reach them, at the 
height of minimum rise. 

For the purpose of calculating minimum capil- 
lary rise® the sand is treated statistically and the 
soil grains moved so that they are in regular 
array with grain centers spaced at a distance 
(2r+d); r is the grain radius and d a spacing 
constant. The value of (2r+d) is given by 


(2r-+d)* = 4(23)r3/3(1—P). (19) 


The funicular masses of liquid are now in the 
form of single rings of liquid wrapped around the 
axis of each pair of spaced grains as shown in Fig. 
2b. We shall consider a single funicular ring and 
determine its size; from the particular value of 6 
required for coalescence we shall determine 
minimum rise and hence the lower limit of the 
funicular zone. 

If we apply Eq. (2) to the new type of ring we 
find 
Ap2rr sin @{(2r+d)—2r cos 6} 

=2(2zr sin @)o cos 8, 


from which if we set (27+d) /2r=q we find when 
Ap=pgh: 


g sec 0=1+0/(pghr). (20) 


Eq. (20) gives the size of the “‘funicular”’ ring at a 
height above the free liquid surface, where @ is 
its angle defined, as before, for rings in contact. 


§ Smith, Physics 4, 184-193 (1933). 
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From Eq. (19) we find 
gq = (2r+d) /2r=0.9045/(1—P)}!. (21) 


Two rings of the spaced hexagonal soil will 
merge when their axes make an angle of 60°, 
hence @ has a value of 30° (cf. Fig. 4d). From 
Eqs. (20) and (21) we find with @=30° that 
minimum rise is given by 


(pghr/o) =1/{1.044/(1—P)!—1}. (22) 


Values of (pghr/a) which lie in the funicular zone 
are, therefore, such that 


1/{1.044/(1—P)!—1} <(pghr/o) <6.452. (23) 


The value (pghr/a)=6.452 is, for actual ideal 
soils, close to the value of this quantity at the 
height of maximum capillary rise.° 

Whether liquid enters the saturation zone by 
direct rise or by condensation the final result is 
the same. When the funicular rings of the 
hexagonal soil merge, each elementary tetra- 
hedron becomes filled; additional liquid must 
condense until the meniscus has, in each case, a 
curvature required by the vapor pressure of 
minimum capillary rise; or less depending on the 
position of the tetrahedron with respect to the 
free liquid. Coalescence of these rings will cause 
liquid to enter all cell types. Because of the 
symmetrical nature of the piling the whole soil 
body will become saturated up to the height 
where the rings are no longer able to merge. 
The small distortions required to transform the 
hexagonal soil into an ideal soil will not alter this 
behavior to any degree. Hence the least value of 
(pghr/o) at which saturation sets in, in a soil 
packed to porosity P is given by Eq. (22). This 
behavior may be observed experimentally with a 
few metal spheres, for example BB shot, piled in 
regular array; if kerosene is added carefully the 
mechanism of ring coalescence and saturation is 
easily observed. 

The funicular zone, as stated earlier, is a zone 
of hysteresis; a complex of many completely 
saturated pores, as will be shown later, represents 
the maximum amount of liquid retained, while 
the distribution in funicular form (webs and 
rings) represents the minimum. Exact analysis of 


* Cf. Smith, Foote and Busang, Physics 1, 18-26 (1931). 
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the region is difficult; for it requires a solution of 
the packing problem, as well as a calculation of 
the volume included between the bounding 
surfaces of the actual funicular liquid masses. For 
the purposes of this paper, we shall consider the 
average packing and its distribution as sufficient 
for approximate calculation of the minimum 
retention. We shall regard the “‘funicular’’ rings 
of the hexagonal soil as representing, above the 
minimum height of capillary rise, the normal 
distribution of liquid in the funicular region of an 





v/(2nr*) = a*(q—(k+a) sin {g/(1+«)} J—(1/3)(q—1)?(g+2). 


From Fig. 2b we see that 


(1+a)=q sec 8; 


(k+a) =q tan @; 
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actual soil; the amount of retained liquid is a 
minimum. 

To find the volume of the ring in the spaced 
hexagonal packing, we proceed as before, and find 
expressions for w; and w, identical with those 
given in reference 4. We set, as before, a= R/r 
and k=y/r and find an expression for v/(27r*) 
identical with Eq. (5). We set (2r+d)/2r=q and 
find from the geometry of the ring (cf. Fig. 2b) that 


(1+a)?=(k+a)*+@. 
Applying this relation to Eq. (5) we find 


From Eqs. (24) and (25) we find the alternative form, 


v/(2mr*) =q(q sec @—1)?[1—(x/2—8@) tan @]—(1/3)(q—1)*(q+2). 


From Eqs. (20) and (25) we see that 


a=a/pghr; k= 
and we have finally 


v/(2mr*) = (o/pghr)*\q—[20/pghr+ (o/pghr)? 


(24) 

sin-' {g/(1+a)} =(7/2—86). (25) 

(26) 

{20/pghr+(o/pghr)?+1—q?}!—o/pghr (27) 
+1—q¢ ]}' sin“ [g/(1+¢/pghr) }} —(1/3)(g—1)*(q+2), (28) 


for the volume of the ring of liquid at a height 4 above the free liquid in spaced hexagonal packing. 
qg is given in terms of the porosity P by Eq. (21). For g=1, that is close hexagonal packing, Eq. (28) 


reduces to Eq. (10). 


The expression for the average ring volume @ in the funicular zone can be computed as before; we 


find 


B/(2mr*) = (—agq’/pgr) | F(h2) — F(hi)} /(he— hs) — (1/3) (q—1)?(g+2), 


where F(h) is given by Eg. (13). The required 
values of @; and 6, are obtained from the known 
values of h, and he by use of Eq. (20); the value of 
q required for Eqs. (20) and (29) is determined 
from the known value of P by use of Eq. (21). 
To compute the total retention we must know 
the number AN, of funicular rings per unit 
volume. The packing used is an average one and 
we must therefore consider each pair of adjacent 
grains to support a ring. Each single grain is 
surrounded by twelve others; hence there are six 
funicular rings associated with each grain. The 


(29) 





number of grains per cc of packed space is 
(1—P)/(4/3)xr*; hence 

N;=9(1 —P)/2zxr'. (30) 
The retention Q; per cc of packed space is given 


by 


0,=5Ny, (31) 


where N, is given by Eq. (30) and 6 by Eq. (29). 
The total liquid retained in a volume V of the 
funicular region is 2; V. 

To calculate the distribution, we proceed as in 
the case of the pendular region (Eq. (21)) and set 
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Q;=vN,, (32) 


where v is given by Eq. (28) and N; by Eq. (30). 
It is not difficult, if we neglect distances of the 
order of a grain radius, to compute the total 
quantity of liquid in a saturation zone. The total 
quantity of liquid contained in a cubic centimeter 
of soil packing is equal to P cc where P is the 
porosity; that, in a volume V of packing, is VP. 
Our discussion has rested, so far, on the 
supposition that a free liquid surface exists just 
under the soil body; the purpose has been to 
ensure the existence, in the soil under discussion, 
of all three distribution zones. The effect of 
varying the position of the free liquid surface 
with respect to the soil body, the effect of 
complete removal of such surface and finally 
hysteresis phenomena remain to be discussed. 
The position of the free liquid determines the 
vapor pressures and curvatures existing in the 
soil body; secondarily it determines the types of 
distribution zones present. The free liquid surface 
may be so far below the bottom of the soil column 
that only the pendular distribution of single 
rings exists; as it moves closer to the bottom of 
the soil, it reaches a position that permits 
pendular and funicular but no saturation zones to 
exist; and finally it becomes sufficiently close, so 
that all three zones, as discussed previously, 
exist. The length of the soil column is of im- 
portance; it may be so short that it contains only 
a funicular or saturation zone depending, of 
course, on the relative position of the free liquid. 
There may be no free liquid surface present at 
all; the system is supposed closed and contains 
only the liquid and its vapor distributed through- 
out the soil body. We suppose equilibrium exists 
between the liquid and its vapor. To find the 
distribution of liquid it is necessary to know the 
vapor pressure distribution in the column. This is 
fixed, however, if we know the value of the vapor 
pressure at one point of the column. We shall 
discuss the problem for two positions of the 
reference vapor pressure; first, when it is known 
just beneath the soil column; second, when it is 
known just above the soil column. 
If p is the value of the vapor pressure at any 
point A, either interior or exterior, of the column, 
including the reference points selected above, we 
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note it is possible to place a free liquid surface at 
some distance hy below the given point A 
without disturbing the vapor pressure p already 
existing at A. The two systems, soil column and 
free liquid surface, can be united so that the 
equilibrium of the column is not disturbed. The 
distance hy is chosen so that the value of the 
vapor pressure at a distance /) above the free 
liquid surface is p. This device permits us to 
adapt the previous equations without much 
difficulty. 

We suppose, first, the vapor pressure pp, at a 
point B, just beneath the soil column, is known. 
The vapor pressure po just above the free liquid 
surface placed at a distance ho below, is known 
(po depends only on the temperature). We let h’ 
be the distance from B to the given liquid mass, 
for example, a pendular ring. The height / of this 
mass above the free liquid surface is h=h’+No. 
To apply Eqs. (4) to (32) we need only to use 
this new value of h. We replace o/pghr by 
o/pg(h’+ho)r. There is a pendular distribution, 
for example, when pg(h’+/)r/o > 6.452; satura- 
tion when pg(h’+ho)r/o<1/{1.044(1—P)!—1}; 
and funicular distribution for intermediate 
values. We can, to a good approximation, calcu- 
late ho. We let Dz and Dy be the densities of the 
vapor at the point B and at the free liquid sur- 
face, respectively. The height 


hee f ’ dp/(gD»); 
PR 


we assume Boyle’s law and find Dg=Do(pp/po); 
hence 


ho=(po/Dog) loge (po/ Pr). 


The value of h to be used in finding o/pghr or its 
reciprocal, for the case of a vapor pressure known 
at the bottom of the column, is then 


(33) 


h=h'+ho=h'+(po/Dog) loge (po/ps). (34) 


Secondly, we suppose that the vapor pressure 
pr to be known at a point T at the top of the 
column; we let h”’ be the distance down to the 
liquid mass under consideration, and hy the 
distance of the reference point 7 from the free 
liquid surface. We have 


h=hyp—h" = (po/ Dog) loge (po/pr)—h", (35) 
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to be used in calculating o/pghr for Eqs. (4) to 
(32). The limits of the three distributions can be 
determined as before. We are thus enabled to 
compute the distribution of liquid in the closed 
soil column, provided we know the vapor pressure 
at a single point, exterior to the column. And 
likewise if we know the vapor pressure at an 
interior point; the analysis is similar to that just 
given. 

The distribution above any meniscus of capil- 
lary rise must be discussed, for elsewhere® the 
height of rise has been shown to be variable 
depending on how it is produced. If liquid reaches 
equilibrium by falling, the height of rise is a 
maximum. The value of pghr/o just over the 
capillary surface at this height is roughly that 
for the lower limit of pendular distribution; for 
the value of this quantity at maximum rise is, for 
most ideal soils, close to 6.452, the value given 
earlier for the boundary of pendular distribution.°® 
The height of maximum rise is a stable position of 
the capillary meniscus and will remain so, as 
long as the’free liquid surface associated with it is 
undisturbed. The distribution of liquid above it 
is pendular, except perhaps close to the meniscus 
where it may or may not be funicular depending 
on the value of pghr/o. The height of maximum 
rise marks the upper limit of saturation; above it 
the boundary will not sustain a meniscus ex- 
tending continuously across the pore space. If, on 
the other hand, liquid is allowed to rise into the 
soil it comes to equilibrium at a lower height 
above the free liquid; the distribution above the 
height of minimum rise is at first funicular and 
then becomes pendular. Finally the liquid may be 
forced to an intermediate position; the distri- 
bution is, just above this height, generally 
funicular and, higher up, becomes pendular. It is 
obvious that the total quantity of liquid retained 
in the soil is different for these cases. It is greatest 
when a meniscus is in the position of maximum 
rise and least when in that of minimum rise; it is 
assumed, of course, that the position of the free 
liquid is, in all cases, the same. This is one phase 
of hysteresis phenomena, namely saturation 
hysteresis; there are, because of the variable 
position of the meniscus, several distributions 
and values of the total retention possible for the 
same position of the free liquid. 

Saturation hysteresis is dependent on the 
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manner of filling or draining the soil, and can be 
experimentally demonstrated in the U-system of 
drainage shown in Fig. 1b. In this apparatus the 
liquid may be allowed either to rise or else fall to 
equilibrium. Drainage may be performed under a 
constant value of pghr/o equal to that of maxi- 
mum capillary rise. The meniscus may be left in 
the soil or entirely removed; if it is left in the 
soil the distribution is that above the maximum 
rise meniscus, if it is removed the distribution is 
that given for the minimum meniscus. The final 
distribution is different, although the position of 
the free liquid is the same in each case. Fig. 6 
shows a schematic diagram of the retention cycle. 
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Fic. 6. Schematic diagram of retention cycle: ABCDE 
represents the distribution with the meniscus at minimum 
capillary rise; ABFDE represents the distribution with the 
meniscus at maximum rise; the area BCD F is the region of 
saturation hysteresis. 








With the saturation zone ending at the minimum 
height of capillary rise, the curve ABCDE 
represents the distribution; it is the normal one 
and the retention is a minimum. When the 
saturation zone extends to the maximum height 
of rise, ABFDE represents the distribution. The 
arc CD represents the normal funicular distri- 
bution and is computed from Eq. (32); the arc 
DE represents the pendular state and is computed 
from Eq. (18). The point B is fixed by minimum 
rise and the point F by maximum rise. The area 
BCDE defines the region of saturation hysteresis 
and any point within may be regarded as a stable 
retention value in the soil. These conclusions are 
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in agreement with those reached experimentally 
by Haines.'° 

There is asecond type of hysteresis phenomena, 
a microscopic hysteresis. It is confined to con- 
figurations composed of a small number of grains. 
It is best illustrated by considering a random 
tetrahedron formed by the centers of four 
neighboring grains and the possible distribution 
of liquid within the grains defining it. We suppose 
the tetrahedron to be in the funicular region. 
There are, in general, three possibilities. First, we 
may have simple rings around each actual 
contact point (the tetrahedron is not supposed to 
be regular). Second, we consider a single face 
defined by the centers of three adjacent grains; 
liquid may be confined within the pore formed by 
the grains and is bounded by two capillary 
surfaces, symmetrical about the face, one on the 
outer side and one facing the interior of the 
tetrahedron (this is a simple web of liquid). If a 
single web exists in the tetrahedron, rings may 
exist at the remaining contacts; a second web, 
joined to the first may also exist, with the 
remaining contact occupied by a ring. Finally, 
the whole tetrahedron may be filled; each face is 
covered by a capillary surface similar to the 
outer of the two bounding the single web. In any 
of these cases the Kelvin equation may be 
satisfied; the tetrahedron is for each of these 
forms at the same distance, properly chosen, 
above the free liquid and each is in equilibrium 
with the prevailing vapor pressure. The amount 
of liquid retained is, however, in each case 
different; when rings exist, it is a minimum and, 
when the cell is saturated, a maximum. The 
distribution is thus complicated by the presence 
of a number of isolated liquid masses each 
completely saturating the grain configuration 
supporting it. 

There is a definite limit to the size of these 
liquid forms. As the tetrahedron is farther 
removed from the free liquid the capillary 
surfaces bounding the web recede further into the 
supporting cell, until they meet in its plane of 
centers. The liquid mass must then break into 
rings. The extreme case occurs when the cell is 


+0 W. B. Haines, J. Agri. Sci. 20, 97-116 (1930); cf. B. A. 
Keen, The Physical Properties of the Soil, Longmans, 
Green and Co. (1931), chapter III. 
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saturated; each bounding capillary surface may 
recede until it occupies the position of a meniscus 
at maximum capillary rise, the smallest opening 
of the cell, after which, it breaks into rings. There 
is, therefore, just as in the case of saturation 
hysteresis, a definite upper limit to the total 
quantity of liquid stably retained by the soil 
when microscopic hysteresis is present. No simple 
approximation, such as has been used to find the 
minimum distribution (Eqs. (18) and (32)), 
however, appears possible. 

It is not difficult to fix the limits of microscopic 
hysteresis. The extreme configuration is a close 
regular tetrahedron, that is six actual contact 
points. Each contact point is capable of sup- 
porting a ring until these merge when pghr/o 
= 6.452. The tetrahedron may be saturated until 
the bounding capillary surfaces occupy the 
smallest opening in each of its faces, that is, the 
position taken during maximum capillary rise; 
the height at which this occurs can easily be 
established. We suppose the tetrahedron satu- 
rated and in communication with a free liquid by 
a pipe running through one of its grains. It is 
sufficient to consider a single face. We find 
pgh/o=xr/(3'r?—xr?/2); hence (pghr/c)=19.5. 
The value (pghr/c)=19.5 marks the extreme 
limit of any hysteresis and is well within the 
pendular zone. It is obvious that other more open 
configurations are capable of being saturated as 
we proceed downward. In an actual ideal soil, the 
number of close tetrahedrons is not great and 
becomes smaller as the porosity increases. The 
normal pendular distribution should not, there- 
fore, be seriously disturbed in the neighborhood 
of pghr/o =19.5. But as we approach the funicular 
region, and in this region, larger disturbances are 
to be expected. 

The distribution given by Eqs. (18) and (32), 
as well as the distribution zones given in this 
paper are natural ones for absorption. If liquid 
is allowed to condense from its vapor into a dry 
soil the amount of retained liquid is a minimum 
(providing, as is here supposed, surface con- 
densation is negligible) and the distribution is 
that given by these equations. In the drainage 
problem, on the other hand, liquid is deposited by 
the passage of a capillary meniscus through the 
soil, and the amount of liquid so deposited 
depends entirely upon the manner of drainage. 
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After liquid is placed in the soil by drainage a 
temporary adjustment quickly takes place. 


Whatever flow over the grains takes place occurs - 


at this stage; liquid flows to the parts of the 
grains where it can exist stably, for example, the 
simple rings of liquid at contact points. All 
further adjustments are by distillation from 
some liquid masses and condensation on others; 
the final state is controlled entirely by the 
Kelvin equation. It is believed that this latter 
process is responsible for the slow movement and 
migration of soil moisture; rather than any flow 
over the grain surfaces, which according to our 
view occurs only in the initial stage of moisture 
adjustment and then only locally. When liquid is 
placed in a soil by a capillary meniscus moving 
through it, the microscopic forms are in general 
fixed and the only further physical process 
occurring is that of distillation and condensation 
until the curvatures of the minute capillary 
surfaces bounding these forms are adjusted to 
satisfy the Kelvin equation. From what has been 
said, it is apparent that the soil boundary offers 
several types of surface satisfying the governing 
Kelvin equation and just which of these exist in 
the final distribution depends upon the kind of 
liquid bodies initially placed in the soil during 
drainage. These multiple effects are confined, of 
course, to the hysteresis region. In a straight 
absorption process, from vapor, such as would 
take place in an initially dry soil, they do not 
exist and the normal distribution prevails. 
King" determined the distribution and total 
retention of water in a column of uniform sand 
after it had stood, practically closed, for a 
considerable interval. His apparatus, shown in 
Fig. 1c, consisted essentially of a column about 8 
feet in length provided with a small capillary C 
at the top so that atmospheric pressure prevailed. 
The liquid was allowed to drain into the small 
bottle D at the bottom, also provided with a 
capillary E. King filled the column by allowing 
liquid to enter at the bottom and rise to the top, 
and then drained it. He allowed the soil with its 
liquid contents, to stand for a period of ap- 


" F. H. King, 19th Annual Report U. S. G. S., Part II, 
59-294 (1897-1898); for the experiment here cited see 
pp. 85-95. 
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proximately 25 years; drainage into the bottle D 
was recorded continuously during the interval. 
He finally measured the distribution and total 
retention. To determine the distribution, he cut 
the column into sections 3 inches long, and, by 
weighing, determined the moisture content of 
each. His values are recorded as percentages of 
the total dry weight of sand in the section. Data 
are given for five different sands, all graded. 

There are, however, several objections to the 
use of King’s data for comparison purposes. 
Unfortunately he does not state the vapor 
pressures at either the top or bottom of the 
column. Nor does he state if a free liquid exists 
below the soil; water collected in the bottle D 
would serve the latter purpose, but, as to this, he 
gives no definite information. There is a further 
objection to his method of filling the column prior 
to drainage. When liquid is introduced at the 
bottom, as King has done, large masses of air 
are trapped in the soil; normally saturated zones 
will contain, therefore, less than their full quota 
of liquid for a long period of time. Finally, the 
grains are far from being spherical. 

Figs. 7a to 7e, inclusive, show curves for the 
data of King and represent, for each of his 
samples, the percentage of saturation existing 
in a lamina at the height h above the bottom of 
the soil. For these curves his data giving the 
amount of retained liquid as a percentage of the 
total dry weight of a section have been expressed 
as percentages of saturation. It appears reason- 
able, for the purpose of plotting these results, to 
consider the variation of the liquid quantity to be 
linear in each section; hence the percentage of 
saturation for each section is taken to be that at 
its center. Curves showing the distribution 
predicted by Eqs. (18) and (32) are also given; 
these are represented by a full line while a dashed 
line is used to represent the experimental results 
of King. 

It will be observed that each of the calculated 
curves of Fig. 7 shows three well-defined regions: 
a horizontal part AB which represents complete 
saturation; a steep part CD corresponding to the 
funicular distribution (Eq. (32)); and finally a 
flat portion DE for pendular distribution (Eq. 
(18)). The abscissae of the points B and E are 
each equal to the height of minimum capillary 
rise; the point D marks the limit of pendular 
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Fic. 7. Distribution of retained water in King’s experiments. The full lines represent the distribution calculated from 
Eqs. (18) and (32); AB is the saturation region; BCD is the funicular; DE is the pendular; m marks the limit of micro- 
scopic hysteresis. The dashed lines represent King’s data. (a) Sand No. 20; (b) sand No. 40; (c) sand No. 60; (d) sand No. 


80; (e) sand No. 100. 


distribution. BC is a straight line parallel to the 
axis of ordinates and extends over the dis- 
continuity between the saturation and funicular 
regions of the soil. BCD may be considered to 
define the funicular region of the soil. The same 
division is suggested, to a greater or less degree, 
by the experimental curves. 

Equilibrium, it should be observed, is estab- 
lished by distillation of vapor from one mass and 
condensation on a lower one. If the system is 
closed and contains besides the soil only the liquid 
and its vapor, it should be established quickly. 
When, on the other hand, the pores are filled with 
a mixture of air and vapor, the time to reach 


equilibrium is lengthened. The early stages of 
readjustment proceed rapidly, but not so with the 
later stages; equilibrium depends on the diffusion 
of vapor through the air in the pore space, and 
the time to reach the final state may, in packings 
with small grains and hence very narrow chan- 
nels, be, therefore, considerable. Comparison of 
the curves would indicate that the coarser of 
King’s sands are pretty near to equilibrium, 
while the finer are not. Hysteresis, of course, is in 
all probability present, and, though the sands 
were in equilibrium, the experimental curves 
would still lie above those given by Eqs. (18) and 
(32). 
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TABLE I. Summary of values of zone and total retention for King's experiments.* 








Sand number 

Porosity (P) 

Radius of grains (cm) 
Height of column (cm) 


Minimum capillary rise (cm) 
Maximum capillary rise (cm) 
Limiting height of microscopic hysteresis (cm) 


Length of saturation zone (cm) 
Length of funicular zone (cm) 
Length of pendular zone (cm) 


(Average vol. of liquid per cc of saturated zone) /(porosity) 
(Average vol. of liquid per cc of funicular zone) /(porosity 
(Average vol. of liquid per cc of pendular zone) /(porosity) 


V =Total volume of packing (cc) 


V’ = Vol. of liquid required to saturate the total pore space (cc) 


Vol. of liquid in saturation zone expressed as a percent of vol. 
Vol. of liquid in funicular zone expressed asa percent of V’ 
Vol. of liquid in pendular zone expressed asa percent of V’ 


Total vol. of liquid retained expressed asa percent of V’ 
Calc. vol. of liquid retained = 


Weight (grams) of liquid retained (density = 1) 
Obs. weight (grams) of liquid retained (King) 


Calc. weight (grams) of liquid retained expressed as a percent 


of dry weight of sand 


20 40 60 80 100 
0.389 0.401 0.408 0.406 0.398 
0.02373 0.00924 0.00776 0.00592 0.00413 
243.9 243.9 243.9 243.9 243.9 
13.40 33.26 38.80 51.10 74.92 
18.71 45.55 52.29 69.76 102.9 
60.4 155.0 184.6 242.0 347.0 
13.40 33.26 38.80 51.10 74.92 
6.56 18.00 22.70 29.10 39.70 
223.9 192.6 182.4 163.7 129.3 
1.0 1.0 1.0 1.0 1.0 
0.200 0.199 0.198 0.191 0.197 
0.0109 0.0250 0.0285 0.0301 0.0516 
30910 30910 30910 30910 30910 
12030 12400 12610 12550 12300 
V’ 5.50 13.60 15.90 21.00 30.70 
0.54 1.47 1.84 2.28 3.21 
1.01 1.97 2.13 2.02 2.74 
Fol 17.0 19.9 25.3 36.7 
854.0 2110.0 2510.0 3158.0 4530.0 
2121.0 2475.0 3515.0 4576.0 5832.0 
1.71 4.25 5.09 6.53 9.15 
4.24 5.05 7.25 9.41 11.82 


Obs. weight (grams) of liquid retained expressed as a percent 


of dry weight of sand (King) 








* Water is the liquid used in these experiments; its density is assumed, in these calculations, to be unity and its surface 


tension 72 dynes/cm. 


Table I shows the retention in each zone 
together with the total retention, computed by 
the methods of this paper. The average volume of 
water per cc of packed space in the pendular zone 
has been calculated from Eq. (17), and for the 
funicular zone from Eq. (31). All calculations are 
made on the basis of a free liquid surface present 
at the bottom of the soil; the assumption is made 
that this surface exists in the collecting bottle D. 
King’s values for the total retention are also 
shown. 

The experiments of King, it appears, should be 
repeated under better controlled conditions. It is 
advisable to use, for experiments of this kind, 
carefully graded glass pearls such as those used 
by Green and Ampt in their work on capillary 
flow or by Haines.'® A small sample of packed 
soil, capable of being weighed on a chemical 
balance ought to suffice. The container for the 
sample should be so constructed that the soil 
pores can be evacuated and then filled with 


liquid. The system should be placed, during 
drainage, in communication with the liquid vapor 
at known vapor pressure and the moisture 
content determined after equilibrium has taken 
place. The procedure used by Van Bemmelen, 
Anderson” and others in investigations on 
absorption by silica gel would be helpful; for 
example, the device of using mixtures of sulfuric 
acid and water would be useful in securing 
definite vapor pressures, although complications 
might arise since the soil body is initially satu- 
rated. It ought, however, to be possible, with a 
procedure of this kind, definitely to control 
distribution experiments. Finally microscopic 
studies of single pore cells with their associated 
liquid, and made at definitely controlled vapor 
pressures, are desirable. 


2 J. S. Anderson, Zeits. f. physik. Chemie 88, 191 (1914); 
also J. W. McBain, The Sorption of Gases and Vapors by 
Solids, London (1932). 
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